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Аннотация 

Введение. При численном решении задач теории упругости в трехмерной постановке методом конечных 
элементов применяются конечные элементы (КЭ) в форме параллелепипедов, призм и тетраэдров. Обычно 
построение матриц жесткости объемных КЭ базируется на принципе изопараметричности, суть которого 
состоит в использовании для аппроксимации геометрии и перемещений полиномов Лагранжа. В расчетной 
практике наибольшее распространение получили так называемые полилинейные изопараметрические КЭ 
с линейным законом аппроксимации перемещений. Главный недостаток данных элементов кроется в эффекте 
«юосвте» («запирания») при моделировании изгибных деформаций. Причем погрешность численного решения 
существенно возрастает в случае, когда конструкция, по сравнению с обычными деформациями, претерпевает 
значительные смещения как жесткое целое. Многолетний опыт решения задач механики деформируемого 
твердого тела методом конечных элементов показал, что существующие объемные КЭ обладают медленной 
сходимостью при моделировании изгибных деформаций пластин и оболочек. Цель настоящего исследования 
состоит в построении на основе метода двойной аппроксимации матриц жесткости полилинейных объемных 
КЭ повышенной точности, позволяющих учитывать жесткие смещения. 

Материалы и методы. Для построения матриц жесткости объемных КЭ применен математический аппарат 
метода двойной аппроксимации, суть которого состоит в раздельном представлении функций распределения 
перемещений и деформаций внутри элемента. Хранение и обработка результирующей системы уравнений 
реализованы в алгоритмических терминах разреженных матриц. Разработка программного обеспечения 
и проведение вычислительных экспериментов осуществлены с использованием 64-х разрядной вычислительной 
платформы М1сгозой У15иа1 За 2013 и компилятора пЁе® РагаПе| ЗаАю ХЕ 2019 со встроенным текстовым 
редактором Пие!® \У15а1 Еопгап Сотрозег ХЕ 2019. Визуализация результатов расчетов выполнена с помощью 
дескрипторной графики пакета компьютерной математики МаЙаб. В качестве тестового образца использован 
объемный восьмиузловой КЭ $ОТЛО185 программного комплекса АМЗУ$ МесВатса|. 

Результаты исследования. Разработано математическое и программное обеспечение для исследования 
напряженно-деформированного состояния массивных конструкций при различных видах внешнего 
воздействия. На тестовых примерах с известными аналитическими решениями выполнена верификация 
авторизированного пакета прикладных программ. Показано, что построенные КЭ по точности удовлетворяют 
основным требованиям, предъявляемым к конечно-элементному моделированию пространственных задач 
теории упругости. 

Обсуждение и заключение. Проведенное тестирование разработанного математического и программного 
обеспечения показало, что построенные на основе метода двойной аппроксимации конечные элементы успешно 
конкурируют с аналогичными объемными элементами $ОГ/ЛО185 программного комплекса АМЗУ$ Месвашса]. 
Предлагаемые элементы могут быть интегрированы в отечественные импортозамещающие программные 
комплексы, реализующие метод конечных элементов в форме метода перемещений. 


Ключевые слова: метод конечных элементов, моментная схема метода конечных элементов, метод двойной 
аппроксимации, объемные конечные элементы, тестирование конечных элементов 
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Е. 1. Уоатевлс ЕЁ: 
а — ас -поде; Ь — $1х-поде 


ТБе геЙайоп$р Бебмееп е соуапапё сотропеп оЁ Фе зат {епзог ш фе юса| Баз15 ап4 415р!асетеп5 ш Фе 
21оба| ахез Ваз Ше Юг [7]: 
1 
5: == ( 2 т, ии: +2 м 


2 
(зиттайоп оуег гереаёе4 ш4ех) 


УВеге 2„,=02„/Ох,; и„, =ди„/дх,, 5=Ь]. 


МесВап!с$ 


367 


Бйр://уезиК-доп$а.ги 


Сийазпигоу РР, ей 1. АррПсапоп ор йе Роие Арргомтайоп Мешо4 юг Сопягисип? 5 тез; Майтсе$ о} Иштештс Етйе Иетет5 


Тре ге]аноп$юр Бебмееп уесюг оЁ ЧеогтаНоп$ {=} ап уесюг о по4да| Чзр1асетет {и} 18 гергеземе4 шт 


табчх Ююпи: {8} =[О2] {и , 

\Пеге Моск тавчх: [Р] = [2 [21,...[Р1,,. ] 
(6хп,) 

забтах: 


[2], =| {2}, {02°}, {02°}, |, =1,2,...л,. 


(6х3) 


Ехргез$$1юп$ ог йе уесюг со[ати$ ог е сопзАегед ЕЁ Вауе Фе ЮПо\уше Югт [10]: ап ею -поде еетепе 


(Е1ю. Та) 


Ри [ил (но +2 Рид нь +2. Рэ)хХз + 


Е аа +212 Рзк 12 из Ри 121 Рэь Ри) хх 


Рок | 2т,2 Е +22 Ри)х, на +22 Ри)х: + 


Рзк| 2т,3 + (2 таз +2 из Ри Хх! + (2.23 + 2 из Рок )хо + 


а" Е. — и в и Я : 
ОС: (2 илоз 2 аз Рок из Ри +23 Ри Рь) хх, | й 


211 Рэк +2 п,2 Ри + (2 дз Рэк +21 Рэк Рзк 1 2т,23 Ри + 


+2 „Ри Рук): 


211 Вэ + 2из Ри (Ен Рак 2 1 Рок Рук + 2из Ри + 


+2 „з Ри Рок) х› 


2 то Рук Т2тз Рок (нь Ру +22 Рик Рак +2 из Ро + 
+2 „з Ра Рак )х, 


1х-по4е @етепе (Е1е. 1 Ь) 


ил Фил + (из Фа +21 Физ) № 


‚2 Фё2 + (2 из Фе +2 Физ) хз 


2 из Физ на Физ Физ) х, + 


я 1 +(2 „3 Физ +2 „з Физ) > 
{2 + 8 г 


РФ Фы (Ев ФЕ ФВ + 


ТЬе поайоп 1$ шиодисе4 Беге: 
02 Р-Я 


— т — 
Е. = Е 
"“  дх и дх.дх, 


х=Хо=Х. =0 
3 
92 


Е = 
т,123 
дх, дх, дх, 


2 


ж=хХо=хз=0 


ХЕХ =хз=0 


_ ОФ, (я > 9°ф, 
Фра _дх, ) Фев > 


жЕХо=хХз=0 


; о, В =1,2,3. 


дхдх, 


жЖЕх.=хз=0 


(2) 


(3) 


Ехргез1опз Гог “зВаре йапсйопз ф, (х,х,,х,) оРа их-поде ЕЕ оМатеа оп ве Баз15 оР Ве ро[упопиа] (1) изте 


Фе дезепеганоп рипстр!е Бауе Фе ЮПо\уше Ююгт: 


1 
ф, НР Рыхы 


ы 1 
Фф, РО ры Р 


А4уапсей Епэтеетпиз Кезеагсй (Ко51оу-оп-Ооп). 2023;23(4):365-375. е155М 2687—1653 
= 


5 1 
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ХЕ сошр!ег уу Фе БаШ-шт Пие® У15 а] Еопгап Сотрозег ХЕ {ех( е4йог. ТБе ргосеззез о# $юогие ап ргосезз те 
Фе 21оБа| $Ипез$ плаблх \еге паретеще@ ш {егп$ оф зрагзе тайлсез [11]. То у1зча|хе Фе гези!$ оЁ са1сшайопз, 
Фе дезспрюг этарысз оР фе МАТГАВ сотршег зу$ет \аз изеа. 

Везеагсв Кези 5. Те шуезисаноп оп Ше ассигасу ап4 сопуегоепсе оЁ Фе 4еуе!ореа йпие @етепе а1гог та 
уаз саге оп оп {её ехатр[ез \Ий ап апа!уйса| зоайоп. ТВе {е5{ ехатр!ез зВо\ питепса] зоийоп$ оМате4 из 
Фе 4еуеюре@ еетепз ап@ Фе 50110185 @етепё оЁ Фе АМЗУ$ Месрапса| зоЙ\аге расКазе зпаИаг ш 
Аитепз1оп [5, 11]. Ве!о\ аге ехатр!ез з@ес{е4 зо {тай Шеу сощат а сотбтайоп оЁ Бепате 4еогтайоп$ апа пела 
Ч15р!асетеп$ оЁЕЕ. 

Ехатрие 1. А зрИсе4 гте полу Ихе4 ш опе зесНоп апа 1оа4е4 уу сопсепгае4 Югсе аё Фе Нее спа. ТВе 4е$1юп 
зспете оЁ Фе ге 1$ зпомп ш ЕР1еге 2. ша даа: ауегасе гадтл$ А = 0.2 т; сгоз$-зесйопа] Аптеп$1юн$ 0.2х0.2 ст; 
тоди!а$ оЁ е!азисиу Е = 101 Н/п?; Ро15зоп'$ гаНо у = 0.3; сопсепгайе4 Ютсе Ё= 10 Н. 

Тве деЙесноп а Фе рошё оРаррИсаНоп оЁ Рюогсе ассог4ше №0 йе еогу оЁ сигуе4 го@$ (ехас{ зо[аНоп) 1$ [7]: 


Ел? —10-3.14-0.23 
= = 0.00302 т. 


ЕЛ 1-10" .8.333.107° 


Ее. 2. Сасщаноп зсБете оЁ®е зрИсе4 гтя 


ТБе сопуегоепсе гези!65 аге ргезеще ш ТаЫе1. 


ТаЫе 1 
Сопуегоепсе гези$ Рог фе зрИсе4 ге 
ана 1/6, м/% 
РАМ АМЗУЗ 
2х2х32 —0.00170 / 44 —0.000433 / 86 
2х2х64 —0.00292 / 3.3 —0.00125 / 59 
2х2х128 —0.00293 /3 —0.00251/ 17 


ТаЫе 1 зВо\уз Ше уаше ог 4еНесйоп /ш Фе патегаюг, ап4 гейаНуе еггог 6 ш Фе депоплтаюг. 


Те у1зиаН2айоп оЁ 41зр!асетеп и, оМате4 оп Фе Баз1 оГРАМ ап4 АМЗУ$ оп а 2х2х64 си4 1$ зпо\уп Ш 
Е1оигез 3 апа 4. Мо{е Ша Ше Не!4 оЁ уегиса| 4зр]асетеп оатей изше АМЗУЗ 40оез поЁ геНесё Фе 7опез %ИВ 


ехаете уаше и ,=-—0.00314 т зВо\уп ш Е1сиге 3 Бу Ше ато. 
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и, т 
-0.003138 
-0.002456 
„^ -0.001793 
-0.001121 
-0.0004484 
0, 


Ес. 3. О15ёайоп райет и, (РАМ) 


и, м 
-. 001346 
-. 001186 
-. 001026 
м а х м -. 8656-03 
7 -. 705Е-03 
-. 5456-03 
-. 3856-03 
-. 2256-03 
-. 646Е-04 
.956Е-04 


Ев. 4. Р1зблби@оп рацег И, (АМЗУ5) 


> 


№1 115 ехатре, 1е ппропапсе оР{аКтз ш®ю ассоип пела Ч15р1асетепи 1$ рагаса у сеайу зееп. 
Ехатрие 2. А зацаге р1айе, полу ршеБе4 аопэ фе сопютг апа |оа4е ми а ип огоу Фет ще4 юа4. Зоигсе 
даа: з14е 1еп? стороны @= 1 из; 1сКпез$ д = 0.01 ип; е!азисиу тодши$ Е = 10$ Н/т?; Ро13301'5 гано у= 0.25. 


Те ехасё уаше оЁ 1е 4еНесНоп ш Фе сещег ое р|ае 15 даегттеа Бу Юга [12]: 


4 
да 
. 
ЕЙ’ нь в амшы И 
УПеге © = 0.00126; р = РУ) суппанса| Иез; 4 = 0.00888889 НЛи” — а1зтоиеа 1оа4 ищепзНу. ЕхасЕ 
—у 


уаше и, (т теегз) 15 едиа! ю соеЁйслепе 0%. 

ш 115 ехашр!е, № рай оР е ре \уаз сопу4егей факте ш®ю ассоипё е соп@юопз$ оЁ зуштейу. ТВе 
сопуегзепсе гези$ т @е Фогт о этар$ и, -5 юг зтёе-1ауег ап4 дон е-!ауег то4е! аге зпо\уп ш Егез 5 
апа 6. Неге апа Раг®фег, 1ауегз теап Фе бгеаКао\уп ог Фе рае шю ЕЕ Бу 1сКпез$. ш езе Ноигез, уаез 1, 2, 3, 4 
О рагатеег 5 1,2, 3, 4 сотезропа 0 Фе 21145: 4х4, 8х8, 16х16, 32х32. ТВе этарвз Бео\у зВо\ Ше гези5 оЁ Ве 
зошНоп оМатей изшя АМЗУ$З (Ппе 1) апа РАМ (Ппе 2). ТЬе Вопхота1 Ппе ш@саеа Бу Фе питбег 3 соггезроп4$ 


{0 Ше ехас( зоаНоп. 
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5 

-0.0002 -0.0002 

-0.0004 -0.0004 1 

-0.0006 -0.0006 

-0.0008 -0.0008 

0.0010 -0.0010 

0.0012 3 -0.0012 3 2 

и, п Е1о. 5. ОгарВ$ ит 5; Иа, т Е1о. 6. ОгарЮ$ и, 75; 
зшее-ауег зсВете суо-1ауег зсВете 


[ ЮПо\з Нот Е1юиге 5 фаё ИВ а зше-1ауег БгеаКао\уп зсКете, фе гайуе еггог уашез оп а 32х32 ой аге: 
50110185 — 16%; РАМ - 10.5%. \МВеп изше а \о-1ауег зсВете (Е1е. 6) оп а 32х32 виа, уме Бауе: 
50110185 — 36%; РАМ — 2.8%. 

ЕЕ рай {езНпз \/аз регогте4 ог 1е 16х16х2 БгеаКдо\уп зспете \уИй опа @югНоп (Е!. 7). ТВе рак 1е 
гезаН$ т фе огт оЁ4еНесйоп 4156 1Биноп райегиз$ и, аге зно\уй ш Е1згез 8 апа 9. 


Е. 7. ЭсБете оЁ р!айе БгеаКдо\уп г раЬ е5 


и ‚Ш и ‚Ш 

-0.001263 

-. 958Е-04 

-. 851-04 -0.001010 

-. 7458-04 

-. 639Е-04 -0.0007575 

-. 532Е-04 

-. 426Е-04 -0.0005050 

-. 3198-04 

-. 213Е-04 -0.0002525 

-. 106Е-04 

р 0.0 


Ев. 8. О1залбивоп и, (500185) Ев. 9. Рзилбийоп и, (МДА) 


Аз сап Бе зееп Нот Фе Е1еигез, фе ФЧютгНоп оЁ е 214, уВеп изше 50110185, самзез а шоге пойсеае 
азуттену оЁ НеЧ и, , пап ууПеп и515 РАМ. АЕ Фе зате ите, пе уаше оЁ®е тахипит 4ейесйоп ог ЕЕ РАМ и, 
= 0,001263 1 сошс!4ез ИН Ше ехасе зойоп. 

Ехатрие 3. А гоипа рае, ие у ршсВе4 а1оп$ Фе сопютг ап4 |оа4е зу а ап оготу 4154 |оа4. Кад$ 
ап4 ф1сКпез$ оР\е рае: А = 1 т; й= 0.01 т. Те тесБатшса] сопзёап($ аге зпп!аг ю Фе дайа ш Ехашр/е 2. 

ТВе ехас( уаше оЁ Фе деЙесНоп ш фе сещег оЁ Ше р! 15 деегттеа Бу Югиша [12]: 

4 
ве 
* 64 
АЕ Ше уаше ог Фе пиепзИу ог @зёчЬшед 1юа4 4 = 0.00888889 Нлп? уаше и,= 0.01563 т. 


п Фе {е5бпс ргосезз, мо уегз1опз оЁ запрИпе Фе 4 рам оЁ Фе рае (зесюг) оп ЕЁ \уеге изе4. ш Фе Йгз( уегз1оп, фе 
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{гее $14е$ оЁ Фе зесюг \уеге Чту14е шю ап еда! питбег оЁ зеотет. ТВе зесопа уегзоп \’аз Базе оп а гад1а] тееаг 
рамегт оЁ зесюг БгеаК4о\уп. ш 15 сазе, Фе пипабег оЁ е]етеп5 а1опе ве гадта$ ап@ фе сисшаг рай оЁ Фе зесюг \уаз 
аззитеа ю Бе Фе зате. Тве сопз14еге зесюг @зсгейтайоп уегз1опз Юг а 32х32 214 аге зВо\уп ш Е1ютге 10. 


\Уегз1оп 1 \Уегз1оп 2 


ох 
у 
$ А 


Ею. 10. Затр|Цле уегзопз ог фе сисч]аг рае зесог оп ЕЕ: 
1 — ип они бтеаКао\уп оЁ Фе @гее з14ез оЁ Ве зесог; 
2 — гала! тезщаг Бгеак4о\уп 


Те сопуегоепсе гези$ ш фе Югт о эгарй$ и, —5 Юг затрИпз фе % ра оЁ Ше рае ассог4тз 1 уегз1оп 1 


у зшя[е-1ауег ап4 Куо-1ауег БгеаК4о\уп зсВетез аге ргезее4 шт Е1еиге$ 11, 12, апа Е1югез 13, 14, гезресйуеу. 


1 2 3 4 1 2 4 
9 Г . 5 


-0.002 -0.002 
-0.004 -0.004 
-0.006 -0.006 
-0.008 -0.008 
-0.010 -0.010 
-0.012 -0.012 
-0.014 -0.014 
И, т И т 
Е1о. 11. Старз и, — 5 юг уегз1оп 1; Ев. 12. Старз и, — 5 Юг уегзоп 2; 
зше]е-ауег зсВете зше-[ауег зсБете 


ш Фезе Нзигез, уаез 1, 2, 3, 4 о# рагатейег 5 1, 2, 3, 4 сотезроп@ ю Ше ©1145: 4х4, 8х8, 16х16, 32х32. Аз т 
Ехатр!е 1, Нпе 1 сотезроп4$ © а 5ОГЛО185-Базей зоайоп, ап Ппе 2 соггезропа$ ю РАМ. ТВе Воттотиа| Ппе 


шасаеа Бу Ше патЪег 3 соггезроп4$ ю е ехас зо оп. 
[{ ЮПо\з$ Нот Фе абоуе эгарВ$ а ап @етепё сопзгисе4 ассог4ате 0 фе тотепе зсВете оп а 32х32х2 эп оЁ 


уегз1оп 2 Баз а ге]апуе еггог оЁ4 %. 
У1зиаПтаноп рабегиз оЁ фе Не! оЁ @за1ЬиНоп оЁ уегиса| 915р!асетеп5 и, Юг РАМ апа $0110185 \ИВ а 


тада! 1ауоие оРШе 4 рае (32х32х2 вп4) аге зВо\’и ш Е1юиге$ 15 ап4 16, гезреснуе|у. 
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Е1о. 13. ОгарВ$ и, — 5 юг уегз1оп 1; Ев. 14. Старз и, — 5 Юг уегзоп 2; 
с\о-1ауег зсНете с\о-1ауег зсНете 


и, шт и, т 


-0.01503 
-0.01202 
-0.009016 
-0.006011 
-0.003005 
0.0 


Е1е. 15. О1561бивоп и, (МДА) Ес. 16. О1зёиНоп и, (30110185) 


2? 
-. 004536 

-. 004032 

-. 003528 

-. 003024 

-. 00252 

-. 002016 

-. 001512 

-. 001008 

-. 504Е-03 
. 410Е-06 


Егот Фе абоуе Йзигез, # сап Бе зееп фай дезрие Фе диаЩануе сошс14епсе оЁ раЦегиз и. ‚ ге]аНуе егогз Юг фе 


тахнпит еЙесНоп аге: РАМ — 3.8%; $0110185 — 71%. Зисй а $1ет1Йсапе етог, \Упеп изше ЗОГЛО185, 1$ 
Бесаизе Фе 4еуерег$ аррИе “зВаре РапсНоп$” зипПаг ю Фе Ишсйоп$ зе Юг Ше ею -поде ейетеп ю 
арргохппае 1е сеотету ап4 Чзр1асетепт оЁе з1х-по4е е]етепь 1.е., \УИоие Фе “дезепеганоп” рипстре [7]. 

015си51оп апа Сопот. ЗН Ёпез$ тафлсез оЁ уо[итейлс пли пеаг Нпие е!етеп5 сопзгасе4 оп фе Ба$15 оЁ 
Фе а4очЫе арргохипаноп тефо4 таке # роз е ю знашае Фе гез$-5тали зайе о БиПате згасвиге$ оЁ атоИтагу 
зеотену ип4ег уапоиз {урез оЁ ежегпа! асНопз. ТВе Вп4атеща! @1егепсе Бебмееп е ргорозе4 сопсерЁ ап фе 
ргемои$[у Кпо\п ИпИе е]етепе {есбпо]оз1ез 1$ фа йе Ч1зр1асетеп, т 15 сазе, аге её ш э1офа| соот4таез, ап4 
бе сотропеп$ оРе збгат {епзог аге 4еегиите4 т 1оса], ш сепега], попогфогопа] ахез. 

ТБе {е5{ ехатрез зВо\/ фа Фе уолтейлс ИпИе е!етет5, сопзигасеа Бу фе дочЫе арргохипайоп тео4, Вауе 
за е сопуегеепсе ап4 сотрее зиссеззАИШу уу ап е]етеп{ оР а зии!Йаг буре $ОГ. 1185 оР Ше АМЗУ$ Месвашса1 
сотрийпе сотр/ех. 

ТЬе деуеоре@ тафетаНса! зоЁ\аге сап Бе шнодисе шо дотезис ипрой-забз анте зоЙ\аге сотр!ехе$ 
пиретепип? Фе Япие @етепё тео ш фе Ююгт оЁ Фе @5р!асетепЕ тефоч. 
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